Coherent versus sequential electron tunneling in quantum dots 
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Manifestations of quantum coherence in the electronic conductance through nearly closed quantum 
dots in the Coulomb blockade regime are addressed. We show that quantum coherent tunneling 
processes explain some puzzling statistical features of the conductance peak-heights observed in 
recent experiments at low temperatures. We employ the constant interaction model and the random 
matrix theory to model the quantum dot electronic interactions and its single-particle statistical 
fluctuations, taking full account of the finite decay width of the quantum dot levels. 
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Recent experimental studies of electronic transport 
through nearly isolated quantum dots assess the im- 
portance of quantum coherence and the nature of dephas- 
ing mechanisms in finite interacting electronic systems. 
Of particular interest is the Coulomb blockade regime, 
where the thermal energy ksT is much smaller than the 
charging energy Eq necessary to add an electron to the 
quantum dot. In this regime the conductance depends 
primarily on the quantum properties of the dot, such 
as its resonance levels and the corresponding line widths 
due to the coupling between the dot and leads. Electrons 
are allowed to tunnel through the quantum dot whenever 
the charging energy is compensated by an external poten- 
tial and the dot energy levels are in resonance with the 
chemical potential at the leads (small bias limit). The 
tunneling condition can be attained, for instance, by a 
tunable gate voltage V g . In a typical experiment V g is 
varied to obtain the conductance spectrum, a sequence 
of sharp (Coulomb blockade) peaks. 

Sequential tunnelingis the key hypothesis for the stan- 
dard rate equations |3] used to explain the transmis- 
sion spectrum of quantum dots in the Coulomb block- 
ade regime 0, 0- This probabilistic picture neglects 
non-resonant quantum virtual processes, under the as- 
sumption that the resonant decay widths T are much 
smaller than both k^T and the energy separation be- 
tween the quantum dot resonances Se, namely, r <C ksT 
and r <C Se, a condition often met by experiments in 
nearly isolated quantum dots. 

The early experimental data taken from ballistic 
chaotic quantum dots were successfully confronted with 
the sequential theory by using the random matrix the- 
ory (RMT) to model the dot statistical single-particle 
properties 0, Q . More recently, the analysis of the mea- 
sured conductance peak-heights in the Coulomb blockade 
regimeBl Q show significant deviations from this theory 
jH 0, U> indicating that some physics is missing. The 
inclusion of inelastic scattering processes 0, 0, fill 113. 
spin-orbit coupling and exchange interaction |l4Lll5| 
into the sequential approach expand in interesting ways 



the considered physical processes, adding new parameters 
to the description. Unfortunately, these studies achieved 
only a limited success in reconciling theory with experi- 
ment. 

In this Letter we show that quantum coherence, so 
far overlooked, leads to important corrections to the 
sequential tunneling picture 0] and explains some of 
thepuzzles pointed out by the conductance experiments 
0, The importance of coherent processes is justi- 
fied by noticing that while the sequential theory requires 
r <C fcsT, Se, the experiments satisfy those conditions 
only in average, namely, (r) < A = (Se) and (r) < fc^T. 
Since both the decay width T and the resonance spacings 
Se fluctuate, conductance peaks where T is larger than 
kfjT and comparable to Se are not exceptional. More 
importantly, the study of fully coherent transport, as op- 
posed to the sequential tunneling limit, provides a better 
framework to understand the interplay between coher- 
ence and interactions. 

We describe a quantum dot coupled to external leads 
by the Hamiltonian 
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dot 
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We write the chaotic quantum dot Hamiltonian iJdot as 
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Hdot = ~ er,V g )d+d 3 +^Ln(n- l) 
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where d~J creates an electron in the jth eigenstate with 

energy Ej of the closed dot, N = Y^,j djdj is the elec- 
tron number operator in the dot, rjV g is the electrostatic 
energy due to the external gate (as usual, V g is the gate 
voltage and rj depends on the system specifics), and C 
is the effective dot capacitance. Equation J5J) is the con- 
stant interaction model. In chaotic quantum dots ground 
state energy fluctuations due to interaction effects are 
very small in the large N limit 5']. We also do not ac- 
count for spin and exchange interaction, which were re- 
cently addressed in the master equation framework by 
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Refs. [HE3. The electrons in the leads are treated 
non-interacting, namely 



-Pleads — ^ £ fc,a c fc, a C fc,a > 
fc,aSL,R 



(3) 



where c kli creates an electron at the state of wave vector 
k = (2m*£k) 1 / 2 /h at channel "a" either in the left (L) or 
in the right (R) lead. The dot-lead coupling term is 



coupling 



E E 

fc,a£L,R j 



V {k,a),j C k, a d j 



h.c. 



(4) 



The magnitude of the coupling matrix elements V, k a > j 
determine through a Fermi golden rule 0] the electron 
decay width T, or the tunneling rate T/h in the master 
equation framework. For quantum dots in the Coulomb 
blockade regime (r) is much smaller than the dot mean 
level spacing A. 

The conductance through the quantum dot is ex- 
pressed in terms of the interacting system retarded 
Green's function, (i) = -(i/h)Q(t) ({«*,(«), 4(0)}). 
The evaluation of Gf^ (t) follows the treatment presented 
by Baltin and collaborators 0] and generalizes their re- 
sult to cases where the condition r <C Se is not met. 

The retarded Green's function is written as a sum over 
terms containing different (and fixed) number of electrons 
in the dot 



G 5w = -i w E p * <{<**(*)> 4(o)» 



TV 



(5) 



JV=0 



where Pjy is the thermal probability to find N electrons 
in the dot. This probability considers the full set of occu- 
pation numbers {n^} of the -f/dot eigenstates. Equation 
(J5J can be formally solved by the method of equation of 
motion. In practice, the equations do not close unless we 
assume that the number of electrons in the dot does not 
fluctuate, which means that we replace N by its expecta- 
tion value N [J^. This simplification is entirely justified 
in the cases of interest, where e 2 /C 3> max (r, fc^T). 

The matrix representation of the retarded Green's 
function is then casted as 



oc 

g *=eM 



JV=0 



el - H W 

tl -"dot 



£*(e)J (I - n N ) + 

? (e) njvj 
where the quantum dot matrix elements are 
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dot 



= (E J -efiV g + UN)6 iJ , 
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and U is the quantum dot charging energy, namely, 
U = e 2 /C. In Eq. © we define [njv] i3 - = (n») N 5 U j 
as the diagonal matrix whose entries are the canonical 



occupation numbers of the (closed) dot eigenstates. The 
retarded self-energy matrix elements, due to the coupling 
to the leads, become 



[E*(e 



' 1,3 



£ -. 

k,a.eh,R 



•»,(fc,a) *(fe,a)J 



i0+ - e fc . 
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The coupling matrix elements V< k a > j vary in the energy 
scale of £k and hence are practically constant in energy 
windows comprising several single-particle states. We ne- 
glect such variations to write 



where £ fe V i (k a) V, k a)j /(e + 10+ - e fc , a 



(9) 



The energy dependence due to the principal value inte- 
gral is also negligible in the Coulomb blockade regime, 
since there are no open transmitting channels. 
The linear-response conductance is |17| 

G=^g with g = Jde(-^jT R , L (e), (10) 

where /„ is the Fermi distribution function in the leads 
with chemical potential fi. Tr^l is the system transmit- 
tance that can be directly computed from the retarded 
Green's function 



(11) 



Equivalently, the above expression can also be casted in 
the well-known form T R>L = tr {T R G R T L G A ) [T^. 

To this point our approach is quite general. The 
only important approximation we make requires e 2 /C 3> 
max(r,fcsT). Albeit restrictive, the approximation is 
compatible with the Coulomb blockade experiments we 
are interested in. Our approach is reduced to the sequen- 
tial tunneling one 0] in the limit of T <C min (kgT, Se). 
The main improvement is that we naturally account for 
quantum virtual tunneling processes. Those are signifi- 
cant whenever ksT becomes comparable with T, a con- 
dition often met by experiments. Furthermore, both the 
single-particle level spacings Se and the decay widths T 
fluctuate. Even if in average A ^> (r), situations where 
Se is comparable to T are not infrequent. In these cases 
quantum corrections are important. When the condition 
T/Se <C 1 is always satisfied and not only in average, 
corrections to the conductance become indeed negligible. 
This was the limit analyzed in Ref. 0] for the phase 
lapse problem. Note also the contrast with the case of 
elastic cotunneling at the conductance valleys. There, 
the contribution of the off-resonant levels is of order T/U, 
whereas here their contribution is of order T /Se. 

We switch now to the statistical study of the dimen- 
sionless conductance peak heights g max . This analysis 



allows for a comparison between the results of our ap- 
proach, experiments and the sequential tunneling theory. 
The statistical ansatz is to assume that the underlying 
electronic dynamics in the quantum dot is very com- 
plex and hence the fluctuation properties of its single- 
particle eigenenergies and eigenfunctions coincide with 
those of an ensemble of random matrices 0, |]| . Accord- 
ingly, the single-particle levels display universal fluctu- 
ations and their spacings Se follow the Wigner-Dyson 
distribution. Likewise, the decay widths T are Porter- 
Thomas distributed. The inputs of the statistical the- 
ory are the mean level spacing A and the average decay 
width (r). We consider the dot both in the absence of 
a magnetic field (orthogonal ensemble, (3—1) and in 
the presence of a magnetic field B that breaks the time- 
reversal-symmetry (unitary ensemble, [3 — 2). The later 
is the relevant one for comparison with available experi- 
mental data. 

The numerical implementation is straightforward, but 
costly since Eq. © requires matrix inversions for each 
realization. The canonical thermal quantities P/v and 
(rii) N are computed using the quadrature formula ex- 
plained in Ref. [13], already used for quantum dots 
For ksT < A good individual peak height accuracy re- 
quires taking into account at least 30 levels around the 
resonant one. Between 5 x 10 4 and 1 x 10 5 realizations 
were used for the ensemble averaging. The charging en- 
ergy U is taken to be 50A (the results are quite insensitive 
to U, provided U > A). 

The data of Ref. [l| show that at very low temper- 
atures, fc^T A, the conductance peak-height distri- 
bution does not follow the standard random matrix the- 
ory 0. By accounting for quantum coherent tunneling 
we obtain a very nice agreement with the experimental 
distributions. This is illustrated in Fig. ^ for B ^ 
((3 — 2). In the inset we present our results for the dis- 
tribution of 5 max for B = ((3 = 1). In Fig. □ the 
dimensionless conductance peak heights g max are scaled 
to unit mean. We show the peak heights distribution for 
k B T = 0.1A, (r) = 0.1A (solid line) and (r) = 0.2A 
(dashed line). The histogram corresponds to the experi- 
mental result of Ref. Q available only for B ^ ((3 = 2). 
Different dots have different (r)/A, a ratio that can be 
determined from the experimental <7 max . (r)/A ~ 0.1 is 
representative of the analyzed experiments. We find that 
as the ratio (r) /A is increased, the probability to obtain 
small conductances is suppressed in comparison with the 
standard sequential theory (dotted line). This can be un- 
derstood as follows: If a given resonance has small tun- 
neling rates, the contributions due to virtual processes 
through its neighbors will reduce the chance to obtain a 
very small peak. Thus, we expect P(g max = 0) = 0. 

In the early experiment by Chang et al. special care 
was taken to discard from the statistical sample conduc- 
tance peak-heights that did not fulfill V <C ksT. Hence, 
corrections due to the finite ratio T/A are practically 




FIG. 1: Peak height probability distribution P(g max ) for 
k B T = 0.1A and B / ((3 = 2). The same for B = {(3 = 1) 
in the inset. Our theory for (r)/A = 0.1 (solid line) and 0.2 
(dashed line) is compared with the standard sequential tun- 
neling result (dotted line), and the experimental distribution 
(histogram) [J. 



negligible. This might explain why a good agreement 
with the standard sequential theory was found there • 
Note also that as ksT becomes comparable with (r) the 
assessment of the quantum dot temperature through the 
widths of the Coulomb-blockade peaks becomes unreli- 
able, due to the non-negligible T. 

The experimental results of Ref. 1] show another strik- 
ing and unexplained discrepancy with respect to the stan- 
dard rate equations. This is best quantified by the ratio 
between the standard deviation 5g max and the mean con- 
ductance peak heights (g max ), namely 



5g™* _ \/{(g max ) 2 ) - (g max Y 



(9™) 



(f ax ) 



(12) 



In the experiments 8g max is significantly smaller than 
predicted by the rate equations plus RMT. Recent works 
I! El E3 discuss if such deviations can be attributed 
to inelastic processes |2l|. Our approach explains the 
experimental findings in the low temperature regime 
ksT/A <C 1, where inelastic processes are hard to jus- 
tify. In Fig. [21 we show a g for B ^ ((3 = 2) as a function 
of the thermal energy for different values of (r)/A. The 
inset shows a g for the case when B — ((3 — 1). The 
standard sequential theory results Q are illustrated by 
the dotted lines. 

At low temperatures and as (r) /A is increased, our 
o g is significantly reduced with respect to the standard 
sequential theory prediction. For higher temperatures, 
kBT > 0.5A, we obtain larger a g than the measured 
ones. Furthermore, as the temperature increases our a g 
approaches the standard theory result. Similar behavior 
was also recently found by including the exchange term 
in Hdot flil Il5j. However, at high temperatures we ex- 
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FIG. 2: Normalized peak heights distribution width a g for 
S / (I (the B = case is shown in the inset) as a function 
of fcsT/A, for (r)/A = 0.05,0.1,0.2 (dashed-dot, solid and 
dashed lines respectively). Symbols correspond to the exper- 
imental results of Ref. Q for different dots and the dotted 
lines to the standard sequential theory. 



pect a reduction of the peak heights fluctuations due to 
inelasticity and decoherence. 

The suppression of the weak localization peak was re- 
cently used to determine the dephasing time in open 
quantum dots [UEli This inspired Folk et al. to experi- 
mentally investigate the change in the conductance peak- 
height upon breaking the time-reversal symmetry of the 
quantum dots by applying a magnetic field B, namely 



B^0 



(g r - 



'B=0 



(f"> 



(13) 



B#0 



At zero temperature the sequential tunneling theory gives 
a constant a — 1/4. Inclusion of temperature corrections 
and spectral fluctuations give small changes, essentially 
keeping a ~ 1/4 [H03- In Fi S- El 

we show a as a func- 
tion of temperature for different values of (r) / A. Our 
simulations show that a is larger than 1/4 at low temper- 
atures and decreases with increasing fc^T. This behavior 
suggests that a finite ratio (r) /A enhances more effec- 
tively the conductance in the unitary case than in the 
orthogonal case. Since a is very sensitive to the ratio 
(r) /A, particular care must be exercised when compar- 
ing data corresponding to different quantum dots. As in 
the analysis of a g our results suggest that an additional 
physical process is needed to explain the experimental 
data for k B T > A. 

In summary, we have investigated the effect of quan- 
tum coherent processes on the statistics of the conduc- 
tance peak heights. We found that at very low temper- 
atures this leads to significant corrections to the distri- 
bution of conductance peak heights obtained using the 
standard sequential theory. The relevant parameter for 
these corrections is (r) /fesT. Our study also indicates 
that estimates of the inelastic scattering rates and the 
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FIG. 3: Normalized change in the average conductance a as 
a function of temperature for different (r)/A. 



strength of the effective exchange interaction in quantum 
dots using the peak height distributions need to account 
for coherent tunneling in order to be quantitative. 
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